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Aims. We investigate the confinement and long-term dynamics of the magnetised solar tachocline. 

Methods. Starting from first principles, we derive the values of turbulent transport coefficients in the magnetised solar tachocline and then 
explore the implications for the confinement and long-term dynamics of the tachocline. 
, Results. For reasonable parameter values, the turbulent eddy viscosity is found to be negative, with turbulence enhancing the radial shear in 

■ the tachocline. Both magnetic diffusivity and thermal diffusivity are severely quenched, with the values much smaller than the magnitude of 
7-H ' the eddy viscosity. The effect of the meridional circulation on momentum transport via the hyperviscosity becomes important when the radial 

, shear becomes large (larger than the presently inferred value) due to the negative viscosity. The results imply that the tachocline develops too 

■ strong radial shear to be a stationary Hartmann layer. In the limit of a strong radiative damping where the turbulence is active on very small 
| scales (< 10~ 4 Rq), the eddy viscosity can become positive although its effect is likely to be dominated by the hyperviscosity. In comparison 

£s ■ with the momentum transport, the transport of magnetic field, heat, and passive particles is more severely quenched. The results imply that the 
I ' tachocline of thickness is larger than O(10~ 2 Rq), independent of the strength of magnetic fields. In addition, the momentum transport is much 
O | more efficient than the particle mixing in the tachocline, consistent with the observations. 
P ■ 

a ■ 
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1. Introduction 

A consistent theory of transport in the solar interior (in particular, the tachocline) is essential to the understanding of the evolution 
of solar rotation and magnetic fields and the distribution of chemical species. While turbulence is assumed to be absent in the 
interior in the standard solar model (Stix 1989), observations (e.g., Pinsonneault et al. 1989) and numerical simulations (e.g., 
Riidiger & Kitchatinov 1996) suggest that transport in this region — although not so fast as turbulent transport (e.g., such as in 
the convection zone) — should be faster than molecular processes to be consistent with the current rotational profile and surface 
depletion of light elements. Such a modest transport could be due to waves via dissipative processes (e.g., radiative damping 
of gravity waves). Another interesting possibility, which has not received much attention, is that turbulence is present in the 
interior due to a variety of instabilities (e.g., see Spruit 1999 and references therein), but that the overall transport due to this 
turbulence is considerably reduced as a result of turbulence regulation. Our previous works (Kim 2005; Kim & Leprovost 2006) 
have shown that stable stratification as well as shearing by the radial differential rotation in the tachocline can precisely do this 
as the excitation of gravity waves reduces the stochasticity in turbulent flow while shearing enhances the overall dissipation (see 
also Kim 2004). 

The turbulent transport reduction can also be caused by magnetic fields (e.g., Cattaneo & Vainshtein 1991; Gruzinov & 
Diamond 1994; Kim & Dubrulle 2001; Kim 2006). In the tachocline, a strong toroidal magnetic field of the strength 10 4 ~ 10 5 
G is believed to be present, which can easily be generated when a weak poloidal magnetic field is sheared by differential rotation 
in the tachocline. A poloidal magnetic field here could be either of primordial origin evolving on a long evolutionary time scale 
(i.e., slow tachocline), or generated by dynamo process operating on fast time scale of the solar cycle (i.e. fast tachocline) (see 
e.g., Gilman 2000; Petrovay 2003). Thus, magnetic fields can potentially play a crucial role in the transport of momentum, 
chemical species, and magnetic flux on long and/or short time scales. In particular, on an evolutionary time scale, the tachocline 
may be considered as a boundary layer between the uniformly rotating radiative interior and differentially rotating convection 
zone with latitudinal variation (Riidiger & Kitchatinov 1996; Gough & Mclntyre 1998; MacGregor & Charbonneau 1999), and 
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the dynamics of this boundary layer crucially depends on the values of the effective magnetic diffusivity, eddy viscosity, etc. 
The understanding of this boundary layer thus requires the prediction of these turbulent transport coefficients, derived from first 
principles. Furthermore, the momentum transport and chemical mixing across the tachocline play a crucial role in the evolution 
of solar differential rotation and the distribution of chemical species. In particular, the present solar rotational profile and surface 
depletion of light elements (lithium) (Schatzman 1993) indicate that the angular momentum transport must have been more 
efficient than the particle mixing in the solar interior (e.g. Pinsonneault et al. 1989). This should be explained by a consistent 
theory of the momentum transport and chemical mixing in the tachocline, rather than invoking a crude parameterization as has 
often been done by previous authors. 

The purpose of this paper is to provide a consistent theory of turbulent transport in the magnetised tachocline and then 
investigate its implications for a long-term dynamics of the (slow) tachocline. Special attention is paid to the elucidation of 
different effects on transport of shearing, stable stratification, and magnetic fields, identifying what is most likely to be the main 
mechanism for turbulence regulation in the tachocline. The remainder of the paper is organized as follows. We elucidate the 
effects of gravity-Alfven waves on turbulent transport in Sect. 2. In Sect. 3, we incorporate the effect of shear flow given by the 
radial differential rotation and provide the theoretical predictions for turbulent coefficients in the stratified magnetised tachocline. 
We elaborate on implications of the results for a long-term dynamics of the tachocline in Sect. 4. Section 5 is devoted to the 
discussion of the limit of strong radiative damping. Section 6 contains the conclusions and discussions. 



2. Turbulent transport in 2D MHD with density stratification 

We envision that turbulence in the tachocline is driven externally, for instance, when plumes penetrating from the convection zone 
randomly stir/perturb the region. In the presence of density stratification and magnetic fields, this random stirring will excite both 
gravity and Alfven waves in the tachocline, which will in turn increase the memory of fluid motion which would otherwise be 
random and incoherent. As a result, turbulent transport is expected to be reduced (e.g. see Kim & Leprovost 2006). In Kim & 
Leprovost (2006), we have shown that the stable stratification (gravity waves) leads to transport property in the three dimensional 
hydrodynamic (3D HD) turbulence very similar to that in the two dimensional hydrodynamic (2D HD) turbulence (without the 
latitudinal dependences) with a negative viscosity in both 2D and 3D. We recall that without stratification the eddy viscosity is 
negative in 2D HD case while it is positive in 3D HD case. Magnetic fields also tend to make the property of 2D turbulence more 
like that of 3D turbulence, leading to a positive viscosity even in 2D (Kim & Dubrulle 2001). Therefore, we model our stratified 
magnetized tachocline by a 2D incompressible fluid in a local cartesian coordinates x and y. Here, x and y represent radial and 
azimuthal directions, respectively. We consider a uniform toroidal magnetic field B = Boy = (0, -d x ao, 0) and represent the 
differential rotation by a large-scale shear flow Uo = Uq(x)$ in parallel to the toroidal magnetic field. Note that a similar cartesian 
2D model for the tachocline was adopted in Kim & MacGregor (2001). For simplicity, we adopt the Boussinesq approximation 
to capture the effect of density stratification and the quasi-linear analysis by assuming that total mass density p = po + p', particle 
density for chemical species n — no + n', vorticity u — a>' = a>' — d x v y - d y v x , and magnetic potential a = do + a' consist of mean 
and fluctuating components, denoted by subscript '0' and prime, respectively. The governing equations for the fluctuations a/, 
a', and p\ - p'/p (p is the mean constant mass density) can be written as follows: 

(d, + Uody)u>' = gd y pi - B d y V 2 a' + vVV + / , (1) 
(d, + U d y )a' = -v x d x a + r)V 2 a' , (2) 

(3, + Uod^pi = —v x +fxV 2 pi. (3) 
8 

Here, v, 77 and p. are molecular viscosity, Ohmic diffusivity and thermal diffusivity, respectively; / in Eq. ([]]) is the small-scale 
forcing driving turbulence; N = ^-g(d x po + pg/c 2 )/p is the Brunt-Vaisala frequency; c s is the sound speed; po = po(x) and p 
are the mean background and constant mass densities, respectively. Note that the typical values of v, p, 77, and N in the tachocline 
are 10 2 cm 2 s -1 , 10 7 cmV, 10 4 cm 2 s _1 , and 10~ 3 s _1 , respectively. 

To elucidate the role of magnetic fields and stable stratification in transport, it is illuminating to examine their effect on the 
diffusion of magnetic flux. To this end, we ignore the large-scale shear flow Uo and forcing / and recast Eqs. ([T])-© as follows: 

d,6j = igkypi + ik y Bok 2 a - vk 2 co , (4) 
d,a = iBokyColk 2 — r\k 2 a , (5) 
N 2 ikyil) , 

d t p 1 = pr - M Pi • (6) 

8 k 2 

Here, Co, a and p\ are the Fourier transform of at', a' and pi. Equations (O and (O give us d r pi/d r a ~ pi /a N 2 /gBo for small 
dissipation, which can then be used in Eq. (01 to obtain 
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By using Eq. (0, we can compute the magnetic flux as 

T x = {v x a) 

d 2 kdt(d,v x (K)a(-K) + v A -(k)<9 f 2(-k)> 



- T -f[(f)-a(b>)]. (8) 

Here, b = Vxa'2, a = 1 +N 2 /k 2 B 2 ) , and J dt ~ r is used where r is the characteristic time scale of turbulence (see, e.g., Gruzinov 
& Diamond 1994). For stationary fluctuations, the flux r x is related to magnetic energy (b 2 ) as BoT x = rj{b 2 ), thereby giving us 
j] T = F x /B as: 

^r^W (9) 

i;<V 2 ) 

where ifj = t(v 2 )/2 is the turbulent diffusivity in the absence of magnetic field (Bo — > 0) and stratification (N — > 0). 777- in 
Eq. (0 immediately shows that magnetic diffusivity is severely quenched for large R m = r/j/rj as either the magnetic field or 
density stratification becomes strong. This is because the excitation of waves (Alfven or gravity waves) increases the memory of 
turbulent eddies, reducing their stochasticity which is essential for turbulent transport. Note that Eq. (0 recovers the 2D MHD 
result as N — > (e.g. Gruzinov & Diamond 1994). 

To elucidate the effect of radiative damping, it is instructive to consider the limit of strong radiative damping. In that case, 
Eq. © can be approximated as N 2 ik y u>/gk 2 ~ fik 2 pi. A similar analysis then gives us the magnetic diffusivity 

7r-— ^r, (10) 

j;(r 2 ) 

where A — 1 + rk 2 N 2 /fik 4 . The contribution from the gravity waves to the transport reduction (the denominator of rjj) in Eq. (TTOb 
is Tk 2 N 2 /fik 4 , which should be compared with tN 2 jr\k 2 in Eq. {9). The ratio of the two is thus rj//j., which becomes very small 
for strong radiative damping (i.e., for large fi). That is, magnetic diffusivity in Eq. <TT~0b is less reduced by stratification in the case 
of a strong radiative damping, as compared to that in the case of weak damping [Eq. ©]• This is simply because the radiative 
damping weakens buoyancy effect. In the extreme limit of /j. — > 00, the effect of stratification disappears with A — > 1 in Eq. ([Tol l, 
recovering 2D MHD result in the absence of stratification. 



3. Consistent theory with the radial differential rotation 

In Sect. 2, we have shown that both stable stratification and magnetic fields can severely quench the transport of magnetic field. 
In this section, we present a consistent theory of turbulent transport in the magnetised tachocline, by incorporating a background 
shear flow Uoy provided by the radial differential rotation. The latitudinal differential rotation is neglected compared to radial 
differential rotation since it is weaker in the tachocline due to thin tachocline thickness (< 0.03 ~ 0.05 of the solar radius) [see 
Leprovost & Kim (2006) for the dynamics and the effect of latitudinal differential rotation]. This shear flow plays a crucial role 
in regulating turbulence and turbulent transport by shearing [e.g., see Kim (2005)]. For simplicity, we assume a linear shear 
flow Uoy - -x&y, where 31 = \d x Uo\ > is the radial shear (or shearing rate), which is assumed to be positive without loss 
of generality. In order to incorporate the shearing effect non-perturbatively, we employ the special Fourier transform for the 
fluctuating quantities <p': 

cf>'(x, t) = — !— f d 3 k$(k, t) exp {i(k x (t)x + k y y + Lz)} , (11) 
(2?r) J J 

where k x ( f) linearly increases in time due to shearing as 

k x {t) = k x (0) + kyftt . (12) 

Here, = \d x Uo\ is again the shearing rate of the shear flow. In this section, we consider the case where the radiative damping 
rate is small compared to the shearing rate. That is, we consider a strong shear limit where = /ik 2 /^ <sc 1. Here, If = l/k y is 
the characteristic length scale of the forcing. Thus, £, (<k 1) is a small parameter, representing the strong shear limit. For typical 
solar parameters fi ~ 10 7 cm 2 s~' and J[ ~ 3 x 10~ 6 s~', this is a good assumption valid for a broad range of characteristic length 
scales of the forcing l f = l/ky I 10 6 ~ 10 7 cm (~ 10 4 /? Q where ~ 5 x 10 10 cm is the solar radius). Note that the opposite 
limit of the strong radiative damping is considered in Sect. 5. 
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For ^ <k 1, the coupled equations ([B-lO can easily be combined to form the following equation for a = 
a exp [fi(kl/3k y ft + k 2 t)/2] 

d T [(l +T 2 )d T a] + r 2 [(l + T 2 ) + B]a - J^/M, (13) 

to leading order in Here, r = j^/jfe, = k x (0)/k y + ftt, y = \Bok y /ft\, B = tf/Bffi* = N 2 /tu 2 A , co A = \B k y \ is the Alfven 
frequency of the mode k y , and / = /exp [fi(k 3 J3k y ft + k 2 t)/2]. 

The solution to Eq. $T3[ in the limit of strong magnetic field y » 1 and stratification N 2 1$} » 1 can be found as 



"'" ^,rJ T1 [(1+t 2 )(1+t2)(1+t 2 +/?)(! +r2 +y 6)]i/4' 



where 2(t) - J <Jt J l ^*f ■ Note that the strong magnetic field limit y > 1 holds on a broad range of length scales of the forcing 
1/ < Lb = 10 1() cm for Bo = 10 4 ~ 10 5 G in the tachocline. To obtain the turbulent coefficients, we assume that the forcing is 
homogeneous with a short correlation time ly; 

</(ki, r,)/(k 2 , ? 2 )> = T f {2n) 2 6(h - f 2 )5(ki + k 2 )0(k 2 ) . (15) 

A long but straightforward algebra by using Eqs. (TJJ)-(|7]i, (TT4b and (TT~5T > then gives us the eddy viscosity defined by {v x v y - b x b y ) = 
-V T d x U = v T ft (b = Vx a'z): 



- IL. f ' £1 

2ft 1 J (2a 



k m 



) 2 



-i 



i 



i 

ft 1 



r. (i6) 



Here, <^(k) is the power spectrum of the forcing defined in Eq. dT5b ; T ~ f d 2 k<f>(li)/k 2 ~ v 2 /t/ is the strength of the forcing 

with the characteristic velocity Vf and correlation time r/; B = N 2 IB 2 yk 2 = N 2 /a>\, o>a = l^o^l; ft = l^f^ol is me radial shear. 
What is the most remarkable about the eddy viscosity in Eq. ( [ToT l is that it is always negative regardless of the relative strength of 
stratification to magnetic field (i.e., B = N 2 Icj 2 a ). In 2D MHD, the Maxwell stress exactly cancels the Reynolds stress to leading 
order while the incomplete cancellation between the two in the next order gives a small positive eddy viscosity (oc 1 /Bfy (Kim & 
Dubrulle 2001). The cancellation of these leading order contributions in 2D unstratified MHD can easily be checked in Eq. ( TToT l 
by putting B — 0, which gives vt = 0. In contrast, in a stratified medium with B + 0, the leading order cancellation in Eq. ( TToT l 
is not perfect, giving a net negative eddy viscosity. This can be shown to be due to the fact that the exact equipartition between 
the kinetic energy and magnetic energy (for pure Alfven waves) is broken by buoyancy, with larger kinetic energy than magnetic 
energy, driving a negative eddy viscosity. A negative eddy viscosity, in a sharp contrast to a positive eddy viscosity, signifies an 
anti-diffusive momentum transport against the gradient of the shear flow. In other words, the overall momentum transport due to 
turbulence accelerates the mean flow, accentuating its gradient rather than eradicating it. A similar tendency of the acceleration 
of the mean flow in the magnetised tachocline can also be due to the direct momentum deposition of gravity waves via radiative 
damping (Kim & MacGregor 2001). 

For clarity, we examine the behavior of vj as a function of B. For large B (strong stratification/weak magnetic field), vj oc 
[— 1 + \B()k v /N\]/ft 2 , showing that the effect of magnetic fields tends to make eddy viscosity positive. In the opposite limit of 
small B (weak stratification/strong magnetic field), vj oc -N 2 /B^ft 2 , whose absolute magnitude |vr| is small compared to 2D 
HD case. This again reflects the tendency of magnetic fields making the eddy viscosity less negative (i.e. more positive). Note 
that for parameter values typical of the tachocline, the cross-over scale from B > 1 to B < 1 is roughly = 10 7 ~ 10 8 
cm. To recapitulate, the result ( TToT l shows a tendency of a negative eddy viscosity in stratified medium despite the presence 
of strong magnetic fields, suggesting that the turbulent transport in the tachocline would amplify the shear provided by the 
radial differential rotation for reasonable values of parameters. We emphasize that the negative eddy viscosity represents the 
amplification of a large-scale shear flow at the expense of small-scale turbulence. The value of magnetic diffusivity, defined by 
{a'v x ) = —TlrdxCk), depends on whether a = (^/3)(1 +B) 3 ^ 2 is larger or smaller than unity. First, in the case a <K 1, which is valid 
for either B > 1 oi B < 1 on scales L < L M = Blft/jdN 3 ~ 10 8 cm (10 11 cm) for B ~ 10 4 G (10 5 G), we can obtain 



T f r d 2 k 0(k) l 
m = 2Tl J (2af^" G ° 



^ /J l^l ^==r. (l?) 



'0 ^ ' ' y 
1 / x 1 / 3 i 

1 ^/s m 1 
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Here, G = (3/^) I/3 r(l/3)/3; ^ = fik^/M « 1; = Tjtf/fl <s 1; T(x) is the Gamma function; T ~ f d 2 k0(k)/k 2 ~ v 2 f /r f is 
again the strength of the forcing with the characteristic velocity vt and correlation time 77. In comparison, for (B ») a » 1, we 
obtain 



T f r d 2 k W 
m = W 2 } W^ Gx 



1 e l/sfa 



2/3 

where G\ = (3/^) 2 ^ 3 r(2/3)/3. In each case, the density (heat) diffusivity ht ((p'v x ) = -firdxpo) is given by 
fj. 

fx T ~ -77 7 . (19) 
»7 

(see also Kim 2006). 

Equation dTTb shows that the turbulent diffusion of magnetic field (t] T ) can severely be quenched by a strong mean magnetic 
field and stratification, proportional to 1/Bg for/J < 1 and to 1/BqN forB > 1, respectively. It is worth noting that rj T oc 1/B 2 
was observed in numerical simulation (e.g. Cattaneo and Vainshtein 1991) of 2D unstratified MHD turbulence (B = 0). As 
the stratification becomes stronger for a fixed Bq with the further increase in B, the diffusion is now reduced as 1 /N 2 [see Eq. 
1181 11. It is important to note that in all cases, magnetic diffusivity has a much smaller magnitude than the eddy viscosity [in 
Eq. (IT6l l1. with a small value of t]t/\Vt\- Note that a similar tendency was also found in the stably stratified shear turbulence 
without magnetic fields (Kim & Leprovost 2006). Specifically, in the strong magnetic field and weak stratification region with 
a < 1 and B > 1, t]t/\vt\ ~ $ 3 (t]/ 'aO 1 ^ 3 -^ 2 INcja ^ L Furthermore, the heat diffusivity yLj in Eq. ([T9l , although larger than 
Tjj by a factor of fi/rj = 10 3 , is yet much smaller than the magnitude of vj. For instance, again when a < 1 and B > 1, this 
ratio becomes /J-t/Wt] ~ ^ 2 J^ INo>a 1 for typical parameter values. These results have very interesting implications for a 
long-term dynamics of the tachocline, as discussed in Sect. 4. 



4. Implications for a long-term dynamics of the tachocline 

The results of the present and previous paper (Kim & Leprovost 2006; Kim & MacGregor 2001) suggest that the uniform rotation 
in the radiative interior is very unlikely to be explained by hydrodynamical means as the momentum transport in stratified medium 
accelerates the mean flow, sharpening the gradient of radial differential rotation that has been created during the solar spin-down 
(see, however, Charbonnel & Talon 2005). Previous authors (Riidiger & Kitchatinov 1996; Gough & Mclntyre 1998; MacGregor 
& Charbonneau 1999), have however shown that a rather weak poloidal magnetic field in the radiative interior can eliminate the 
differential rotation, thereby leading to a uniform rotation therein. In this case, the tachocline can be envisioned as a boundary 
layer where the generation of the toroidal magnetic field by the shearing of the poloidal magnetic field (due to differential rotation) 
is balanced by the diffusion of the toroidal magnetic field while the dissipation of the radial differential rotation is balanced by 
the azimuthal Lorentz force associated with the large-scale toroidal and poloidal magnetic fields, i.e., Hartmann layer (Riidiger & 
Kitchatinov 1996; Gough & Mclntyre 1998; MacGregor & Charbonneau 1999). By adopting the molecular values for viscosity, 
magnetic diffusivity, and heat diffusivity in the tachocline, these previous authors obtained estimates of the tachocline thickness 
and the strength of the interior poloidal magnetic field. In the case of the tachocline with residual turbulence, for instance, driven 
externally (e.g. plume penetration) or internally (e.g. via instability), the values of turbulent transport coefficients, instead of 
molecular values, should be used in the analysis of the Hartmann layer. 

For the clarity of the discussion, it is worth recalling that the Hartmann layer is based on the configuration where a poloidal 
magnetic field is fully contained in the interior, without penetrating into the convection zone above so that the latitudinal differen- 
tial rotation in the convection zone does not leave its footprint into the radiative interior. By representing the latitudinal coordinate 
by z, and by denoting the poloidal and toroidal magnetic fields by B z and B y = Bq, respectively, the major force balance for the 
toroidal magnetic field in the tachocline B Y and the mean shear flow U y due to the differential rotation can roughly be expressed 
as 

d,U y ~ B z d z B y + v T d xx U y - A T 8 XXXX U y , (20) 
d,B y ~ B z d z U y + T] T d xx B y . (21) 

Here, At ~ yUr(^/^V) 2 (^©/A) 2 is the hyperviscosity due to the meridional circulation (see, e.g., Spiegel & Zahn 1992); Q. ~ 
3 x 10~ 6 s~' and Rq ~ 5 x 10 ll) cm are the average rotation rate and solar radius, respectively; A is a constant of order unity. 
Similarly to a positive eddy viscosity, a (positive) hyperviscosity acts to smooth out the gradient of (or damp) a large-scale shear 
flow. 

While the values of turbulent transport (vj, t]r, and At) for the coupled system (f20l>- (f21~b have conventionally been assumed 
to be positive, the results in Sect. 3 show that the turbulent momentum transport in the sheared stratified turbulent tachocline 
is anti-diffusive with negative viscosity (i.e., vt < 0), accelerating the mean shear flow U y (i.e., radial differential rotation). 
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The negative viscosity will amplifies the shear in the differential rotation and consequently toroidal magnetic field since they 
are coupled via the Lorentz force (via B z ). The crucial question is then how the Hartmann layer is maintained with negative 
viscosity which tends to become unstable due to v T < 0. Obviously, if the magnetic diffusivity rj T is large enough to overcome an 
unstable situation caused by the negative viscosity, the coupled system can find a stable stationary configuration. However, this 
is very unlikely since the magnetic diffusion is pathetically small compared to the magnitude of the eddy viscosity, as discussed 
in Sect. 3. Could the hyperviscosity At ~ (Rq/A) 2 (Q./N) 2 /jt due to the meridional circulation then stabilize the system? To 
answer this question, we note that the contribution from the hyperviscosity to Eq. ( f20b is of order Xjjh A while the contribution 
from the eddy viscosity is of order vt/Ii 2 , where h is of order of the tachocline thickness. Thus, the ratio of the two is roughly 
(fi T /v T )(R Q /h) 2 (Q./N) 2 (e.g. ~ (juk 2 /3l) 2/3 (3l 2 /Nuj A )(Q./N) 2 (Ro/h) 2 for a < 1 and B > 1), which can be shown to be small for 
the parameter values typical of the present solar tachocline. However, it is important to realize that this ratio (oc for a < 1 
and B > 1, for instance) becomes large as the shear J{ increases. Therefore, it is plausible that as the shear is amplified via 
the negative eddy viscosity, the effect of the hyperviscosity becomes important and eventually dominates the negative viscosity, 
possibly stabilizing the system. In order for this to be the case, the shear in the radial differential rotation however has to be 
larger than what is observed today. This can be checked by requiring At/Ii 2 vt > 1. For instance, in the case a < 1 and B > 1, this 
demands M/N > (N / fik 2 ) ll2 (B k y /N) 3/4 (hN / RqQ.) 3/2 , where If = l/k y is the characteristic scale of the forcing. Even if we take 
B ~ 10 4 G and a thin tachocline h ~ 10 ~ 3 Rq, Jl>0AN for k y ~ 1(T 7 cm -1 , and J{ > 0.01N for k y ~ 1(T 6 cm -1 , which seem to 
be rather too large to be reasonable (recall N ~ 10 3 s and J[ ~ 3 x 1CT 6 s are the presently inferred values in the tachocline). 

To make this argument more concrete, it is instructive to examine the behavior of Eqs. (f20b-(f2Tb in more detail. To this end, 
we average them over the space (x) with the approximation d : ~ i/Rq and d x ~ l/h and use the values of vj, t\t and pLj obtained 
for a < 1 in Sect. 3. Here, h is again the tachocline thickness. The resulting envelop equations for J[ ~ U v /Rq and B y , in a 
properly non-dimensonalized form, are as follows: 

d t S = iM + [v(S,M) -pr\{S,Mj\S , 

d,M = iS -Tj(S,M)M . (22) 

Here, S = M/Q, M = B,/Q.R Q , v = or[l - 1/ y/l + l/M 2 ]/\S\ 2 , and rj = 1/ VM 4 + M 2 \S | 2/3 . The parameter a appearing in 
v is a rough measure of \vt/t]t\, which was shown to be very large in Sect. 3 because the magnetic diffusion is much smaller 
than the magnitude of the eddy viscosity. Specifically, a ~ {[ilri) l ^(jMlrik 2 ,) 2 ^N 2 j^ 2 J <9(10 8 ) obtained by using the reasonable 
parameter values in the tachocline J[ ~ 3x 10~ 6 s~ l ,N ~ 3x 10~ 3 s _1 , and ky ~ 1CT 8 cm -1 . The parameter B in Eq. (1221 represents 
a rough measure of the ratio of the effect of hyperviscosity to that of magnetic diffusivity, with the value B ~ (RQD./hN) 2 (fi/r]) ~ 
;<9(10 3 )for/z/fi o > 10 3 . 

We note that the values of rjr, vj and jjj are all proportional to the intensity of the forcing T, which is a free parameter in 
our problem. In our non-dimensionlization, this value is fixed by utilizing the observational evidence that the particle diffusivity 
of light elements (lithium) is about Dj ~ 10 3 cm 2 s~' to be consistent with the present solar surface lithium depletion (e.g., see 
Barnes, Charbonneau & MacGregor 1999). Since r\j ~ (j]/D)Dt where D ~ 20 cm 2 s~' and 77 ~ 10 4 cm 2 s~', we impose the 
condition that 777- ~ 10 5 cm 2 s _1 for the parameter values typical of the present sun (i.e. B v ~ 10 5 G, ~ 3 x 10 -6 s _1 , etc). [Note 
that t]t in Eq. (fT7li-(fT8Tl depends on both J[ and B y , which change in time according to Eq. (1221.1 

The results from the numerical simulation of Eq. ( l22l for a = 10 8 are plotted by the cross symbols in Fig. 1, which shows 
\S\ = for various values of B. The results nicely show that the amplitude of S decreases as B (hyperviscosity) becomes 

large, as expected. However, for small B (hyperviscosity), an amplitude of J[ > 10 3 f2 is too large to be reasonable. For this 
large value of M, which is comparable to the Brunt-Vaisala frequency N, the shear flow is likely to be unstable against a shear 
instability (see, e.g., Drazin 1981). As B increases, the amplitude of J{ decreases and approaches the asymptotic value given by 
the dashed line [a/(2B + 2)] 3 ^ 4 , which is obtained by using |5 1 = \M\. The asymptotic value indicates that J[ ~ 100Q for B = 10 5 
while J{ ~ 3Q for B = 10 7 . By using the definition of B, one can easily show that a reasonable value of shear J[ ~ Q is possible 
only for an extremely thin tachocline with Ii/Rq ~ 10 -5 . Since \S\ = \M\ along the dashed line in Fig. 1, the toroidal magnetic 
field also seems too strong (> 10 5 G) to be stable (e.g. against the magnetic buoyancy instability). To summarize, the results from 
our toy model (l22l clearly demonstrate that for reasonable parameter values, a stable stationary Hartmann layer is very unlikely 
in the tachocline. Thus, a large-scale shear flow is likely to become time-dependent (similarly to the behaviour found in Kim and 
MacGregor 2001), or to develop a secondary instability. 

5. In the limit of the strong radiative damping 

The anti-diffusive momentum transport discussed in the previous sections originates from the stable stratification in the 
tachocline. The buoyancy force can however be weakened by a strong radiative damping, as demonstrated in Sect. 2. Thus, 
it is conceivable that if a radiative damping is large enough to weaken the buoyancy force sufficiently, the momentum transport 
may become diffusive, smoothing out the gradient of the radial differential rotation. In this section, we explore this possibility by 
considering the limit of a strong radiative damping where the density fluctuation p\ in Eq. 0) is stationary with (d, + Uod y )p\ = 0. 
In order for this limit to be valid, the parameter ^ = fik 2 /Jl, which has been assumed to be small in the previous sections, is 
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Fig. 1. The plot of the amplitude of S = Jl/Q. as a function of B (oc hyperviscosity) for a - 10 8 . The cross symbols represent 
the results of the numerical simulation of the toy model Eq. d22l i. The dashed line denotes the asymptotic behavior, given by 
[al(2p + 2)fl\ 



no longer small. We are however still interested in the case where the dissipation rate of magnetic field is small compared to the 
shearing, namely, for g n — 7]k 2 /Jl <K 1. Note that this is a valid assumption for the reasonable values of the characteristic scale 
of the forcing If - l/k y > 10 5 cm in the tachocline. 



By using N 2 v x /g + yN 2 p\ = 0, the coupled equations (Hi-© can easily be combined to yield 



8 r + 



—2 

N 



(1 +T 2 ) 2 



[(1 +T 2 )d T a] +y 2 (l +T 2 )a 



iB 
Ji 2 k 



-fir), 



(23) 



to leading order in <K 1. Here, a = aexp [r}(k?J3k y M + k 2 t)/2]; f = /exp [rj(P x /3k y M + k 2 t)/2]; r = k x /k y = k x (0)/k y + Jit, 
y = \B k y /M\, o>a = \Bok y \, and N = N 2 /Jl(fik 2 ). The solution to Eq. (|23]) can be obtained in the limit £, < 1: 



aw 
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(It 



siny[^(r) - ^(Ti)]/(Ti)e _iv ^M-*^ 
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(24) 



Here, <p(T) = f d^i/fin), ^(t) = 1 - [(1 + t 2 Y 2 -N t(1 + r 2 )" 3 + N (1 + r 2 )- 4 /4]/2y 2 , and m(j) = [tan" 1 r + r/(l + r 2 )]/4. By 
omitting the details of the algebra, we here provide the resulting eddy viscosity vr, magnetic diffusivity r]j and heat diffusivity 
fJ.T- 
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(26) 



(27) 



jy- J \z.nf k 2 

Here, a*^ = \k y Bo\, k y is the typical wavenumber of the forcing, and T ~ r/ J d 2 k(p(k)/k 2 ~ v 2 /Tf is the strength of the forcing 

— 2 — 2 

with a characteristic forcing velocity V/ and correlation time Tf. The function G(./V ) depends on N , taking at most a maximum 

2 2 

value ~ 0.37 for N ~ 4 and vanishing as either N — * or N — > oo. Since G(N ) < 1, the eddy viscosity in Eq. (l25T l is now 
positive. That is, turbulent momentum transport in the limit of a strong radiative damping is diffusive, in contrast to the case of 
weak radiative damping considered in Sec. 3. This can be shown to be related to the equipartition between kinetic and magnetic 
energy to leading order, in contrast to the weakly damped case in Sect. 3 where the excess of kinetic energy over magnetic energy 
causes a negative viscosity. As expected, this is because the buoyancy force is greatly reduced by strong radiative damping. The 
results in Eqs. (T25l)-(|2"7|| show that both vj and r\j are quenched largely by magnetic fields while r\j is also reduced by shear 
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(oc £ 2/3 oc JT 2 t 3 ). In comparison, the heat diffusivity in Eq. (1771 ) exhibits a different behavior, being mainly inhibited by buoyancy 
force. This is because the coupling of the density fluctuation to the fluid is feeble due to strong radiative damping. 

The comparison of Eqs. d25l and ( |26b reveals that the magnetic diffusivity is smaller than eddy viscosity by a factor of 
£^ 3 <k 1 (recall = i]k 2 /Jl is a small parameter characterizing the strong shear limit, where If = l/k y is the characteristic 
scale of the forcing). The ratio t]t/vt « 1 is however larger than the value in the case of weak damping in Sect. 3. This is 
an interesting result since this ratio is crucial in the estimate of the tachocline thickness, which has been set by the molecular 
values with rj/v ~ 10 2 in previous works. We emphasize that our predicted value t/t/vt (<k 1) is much smaller that rj/v based 
on molecular values. Furthermore, the ratio of the effect of the hyperviscosity to that of eddy viscosity on the shear flow is now 
roughly given by (Bo/RQN) 2 (RQ 2 k Y /h) 2 (Q./N) 2 J 1 for reasonable parameter values. That is, the effect of hyperviscosity due to 
meridional circulation is crucial in maintaining the momentum balance in the tachocline. 

Based on these observations, we now seek to obtain the estimate on the thickness of the tachocline and the strength B y and 
B z by requiring d, — and by ignoring vj in Eq. d20l . A simple analysis of Eqs. (l20b- (l2TT) . with the help of Eqs. ( 1261 and d2Tb 
and Uy ~ J{Rq, gives us 



2 

1/3 



h/RQ l_/£l 
R Q k y ~ a\n) 

B z T] T 



By m 2 



C > (28) 

(29) 



where f_ = rjkl/M and It = l/k y is the characteristic length scale of the forcing. In a sharp contrast to the previous works 
which use the molecular values for the viscosity and magnetic diffusivity, the result ( 128b shows that the tachocline thickness is 
independent of the strength of magnetic fields. For instance, for k y = 10~ 7 ~ 1CT 6 cm -1 , Ii/Rq = 0.01 ~ 0.05 (for J[ ~ 3 x 10~ 6 
s _1 , N ~ 10~ 3 s _1 , and A ~ 3). This estimate is much larger than h/RQ ~ 10~ 3 by Gough & Mclntyre (1998), and is mainly due 
to the small magnetic diffusivity with the value j]t//ut ~ £// 3 (N/u>a) 2 <k 1. That is, a more efficient momentum transport than 
magnetic dissipation causes a thicker tachocline. Equation ( 1281 implies that the tachocline thickness decreases as the stratification 
becomes stronger, which could have some implications for other stars. On the other hand, Eq. ( 1291 shows that the ratio of magnetic 
field strengths depends on the tachocline thickness and the turbulent magnetic diffusivity 777-. As 777- is directly proportional to the 
strength of the forcing T , we again utilize the constraint on the turbulent particle diffusivity Dj ~ 10 3 cm 2 s _I in the tachocline 
to be consistent with the surface depletion of lithium on the Sun. Since r\j ~ {rj/D)Dj (see Kim 2006) where D ~ 20 cm 2 s _I 
and 77 ~ 10 4 cm 2 s , we obtain r\j ~ 10 5 cm 2 s _1 . By using this value in Eq. d29l and h/RQ > 0.01 ~ 0.05, we can obtain 
BJBy ~ 10~ 8 ~ 10~ 7 . If we take, B y ~ 10 4 G, B z ~ 10~ 4 ~ 10~ 3 G while for B y ~ 10 s G, B z ~ 10~ 3 ~ 10~ 2 G. These values 
of B z are comparable to the previously estimated value. It is important to note that for r\j ~ 10 5 cm 2 s _1 , Eqs. ( f25l ) and ( |2"7T i give 
vj ~ 10 6 cm 2 s _1 and jij ~ 10 10 cm 2 s (» //). Therefore, our results naturally predict a more efficient momentum transport than 
particle transport (recall Dt ~ 10 3 cm 2 s~') in the tachocline without fine-tuning parameters (see also Kim 2006), as required to 
be consistent with the observations. 



6. Discussion and conclusions 

We have presented a theory of turbulent transport in the sheared stratified magnetised tachocline when turbulence is driven and 
maintained by an external forcing (e.g. due to plumes penetrating from the convection zone). Despite the presence of a strong 
(uniform) magnetic field, the momentum transport is found to be anti-diffusive (vt < 0) for reasonable parameter values in the 
tachocline. The magnetic diffusivity {rjj) and heat diffusivity (fij) are found to be positive, but severely reduced as a result of the 
shear stabilization and transport reduction by Alfven-gravity waves with very small values of t]t/\vt\ (<k 1) and (It/\vt\ !■ 
Since the momentum transport is anti-diffusive, the gradient in the radial differential rotation is amplified. As the shear becomes 
strong, the hyperviscosity due to the meridional circulation, which is always diffusive, can become important, counteracting 
the effect of negative eddy viscosity. However, in order to maintain the radial shear J{ in the differential rotation with a value 
comparable to the average solar rotation Q, as presently inferred, the tachocline as the Hartmann layer has to be very thin with 
the thickness h/RQ ~ 10~ 5 for the toroidal magnetic field B v ~ 10 5 G. Here, Rq is the solar radius. Otherwise, the radial shear 
and/or toroidal magnetic field in the tachocline appears to be too strong to be stable against the shear instability and/or magnetic 
buoyancy instability. Furthermore, since the turbulent coefficients depend nonlinearly on J[ and B v , J[ and B v form a nonlinear 
dynamical system, which can exhibit complex time variation. These results point to the possibility that the tachocline may exhibit 
a much more complex (temporal and/or spatial dependent) dynamics than previously thought even in the slow tachocline scenario. 

In the limit of strong radiative damping where the temperature (density) fluctuation is almost stationary, the turbulent mo- 
mentum transport is found to be diffusive, down the gradient. This requires the turbulence to be on very small scales such that 
= iik 2 /J{ » 1. Here, l/ky = If is the characteristic length scale of the forcing, and 31 ~ 3 x 10~ 6 s _1 is the radial shear. Thus, 
the forcing scale If has to be smaller than 10 6 ~ 10 7 cm, which is about 10 -4 Rq. In this case, the tachocline may be viewed as a 
Hartmann layer with possibly stable configurations of the radial differential rotation and toroidal magnetic field. A simple analysis 
by using the predicted values of vj, jij, and rjj suggested that the effect of the hyperviscosity is likely to dominate over the eddy 
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viscosity. However, since the magnetic diffusion is severely quenched compared with heat diffusivity, the momentum transport is 
effectively much more efficient than the magnetic field diffusion. Consequently, the major force balance in the tachocline leads 
to the estimate of the tachocline thickness ( t 10~ 2 Rq), which is larger than the previously estimated value ~ 10~ 3 Rq by Gough 
& Mclntyre (1998). As rj T and fi T depend on the strength of toroidal magnetic fields and radial shear, the tachocline thickness is 
found to be independent of the strength of (toroidal and/or poloidal) magnetic field (c.f. Gough & Mclntyre 1998). The strength 
of poloidal magnetic field in the radiative interior is estimated to be of order 1CT 4 ~ 1CT 2 G for the value of the toroidal magnetic 
field 10 4 ~ 10 5 G, comparable to the previously estimated value. Importantly, the results provide a natural explanation for a more 
efficient momentum transport than particle mixing in the tachocline. However, a severe reduction in the magnetic field diffusion 
could be problematic for the solar dynamo [e.g., the interface dynamo (Parker 1993)]. 

While the discussion of our results is focused on the applications to the present sun, they might also have interesting im- 
plications for other stars. In particular, our predicted values of turbulent transport coefficients have different dependences on 
the strength of magnetic field, stratification, and radial shear as well as on the molecular values of viscosity, ohmic diffusivity, 
and radiative diffusivity. As the values of these parameters vary from one star to another, it would be of interest to explore the 
implications of these results for other stars. For example, it might be possible to utilize the results to infer the value of the radial 
shear or the strength of magnetic field in other stars, and also to gain some insight into the presence of a tachocline-like shear 
layer in those. Of course, in the case of more massive stars, the rotation rate is much faster the sun, demanding the prediction 
of turbulent transport coefficients by taking into account the effect of average rotation (Leprovost & Kim 2006). Ultimately, it 
will be interesting to investigate a consistent model incorporating the spin-down of the sun (and other stars) due to the angular 
momentum loss. Finally, we note that our analysis, limited to 2D, should be extended to 3D, in particular, to study solar dynamos 
(e.g. the a effect). These issues are currently under investigation and will be addressed in future publications. 
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